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Abstract 

If financial markets displayed the informational efficiency postulated in the efficient markets hypothesis (EMH), ar- 
bitrage operations would be self-extinguishing. The present paper considers arbitrage sequences in foreign exchange 
(FX) markets, in which trading platforms and information are fragmented. In Kozyakin et al. (2010) and Cross 
et al. (2012) it was shown that sequences of triangular arbitrage operations in FX markets containing 4 curren- 
cies and trader-arbitrageurs tend to display periodicity or grow exponentially rather than being self-extinguishing. 
This paper extends the analysis to 5 or higher-order currency worlds. The key findings are that in a 5-currency 
world arbitrage sequences may also follow an exponential law as well as display periodicity, but that in higher- 
order currency worlds a double exponential law may additionally apply. There is an "inheritance of instability" 
in the higher-order currency worlds. Profitable arbitrage operations are thus endemic rather that displaying the 
self-extinguishing properties implied by the EMH. 
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1. Introduction 

Arbitrage operations are profitable if a good or asset can be bought for a lower price than that for which it 
can be sold. Such operations are distinct from speculation in that there is little or no capital or risk exposure. 
Arbitrage profits are made by trading on prices that are already posted, rather than by trading on speculative 
guesses about what prices posted in the future might be. The New Palgmve Dictionary defines arbitrage as "an 
investment strategy that guarantees a positive payoff in some contingency with no possibility of negative payoff 
and with no net investment" , Dybvig and Ross ( 2008 ) . 

In the economics and finance literature the exploitation of profitable arbitrage opportunities is postulated to 
eliminate price discrepancies between goods and assets that are in some sense "identical." This postulate, that 
arbitrage operations are self-extinguishing, implies the "law of one price", this terminology having been coined in 
the purchasing power parity explanation of foreign exchange rates, Cassel (1916). A related no-arbitrage condition 
is that of covered interest parity (CIP), whereby the ratio of forward to spot exchange rates for currency pairs is 
equal to the ratio of the interest rates on comparable assets denominated in the two currencies over the forward 
period in question (Keynes, 1923, p. 13). No-arbitrage condition form the bedrock of mainstream theory in finance, 
as embodied in the Modigliani-Miller theorem regarding corporate capital structure, the Black-Scholes model of 
option pricing and in the arbitrage pricing model of asset prices, Ross (1978). 

The present paper deals with arbitrage processes in foreign exchange (FX) markets. The Bank for International 
Settlements (BIS, 2011, p. 5) provides a useful classification of arbitrage strategies in such markets. Classical 
arbitrage aims to exploit discrepancies between actual exchange rates and no-arbitrage conditions, such as that the 
exchange rates for currency pairs are consistent with their cross-exchange rates. Latency arbitrage aims to exploit 
time lags between trades being initiated and FX price quotes being revised. Liquidity imbalance arbitrage attempts 
to take advantage of order book imbalances between different trading platforms. Complex event arbitrage is geared 
towards properties of FX rates, such as mean-reversion and momentum, or towards systematic patterns in the 
way FX response to new information being released. This paper focuses on triangular arbitrage operations, arising 
from discrepancies between the exchange rates for currency pairs and the cross-exchange rates for the currencies 
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involved. To simplify the analysis we ignore the difference between FX bid and ask prices and do not take account 
of the cost of executing FX transactions. For evidence that there are arbitrage profits to be made, net of bid-ask 
spreads and transactions costs, see Marshall et al. (2007) and Akram et al. (2008). 

If the efficient market hypothesis (EMH), formulated by Samuelson (1965) and Fama (1970), applied to FX 
markets, exchange rates would reflect all the information relevant to their determination and there would be no 
profits to be had from FX arbitrage trades. This would involve an "arbitrage paradox" (Grossman and Stiglitz, 
1976): if arbitrage opportunities did not exist, there would be no incentive for FX market participants to monitor 
exchange rates, so profitable arbitrage opportunities could well arise. There is also the "no trade" problem. If FX 
markets were informationally efficient, as the EMH postulates, you would expect to observe periods in which FX 
transactions volumes dwindled when no new information was being generated. Instead FX transactions volumes 
remain substantial, even during quiescent market conditions. 

The microstructure of FX markets is fragmented, implying that information about profitable arbitrage oppor- 
tunities is also likely to be fragmented. The contrast is between the "lit" areas where information on some of 
the electronic trading platforms is publicly available, and the "dark pools", where information on FX transac- 
tions conducted directly between banks and their end-user clients is initially private. BIS data for 2010 document 
this FX market fragmentation. Electronic booking systems (EBS), such as the London-based EBS and Reuters, 
took 18.8% of global turnover; multi-bank electronic communication networks (ECNs), such as the US-based Cur- 
renex, Hotspot FX and FXall, accounted for 11.1% of global turnover; and single-bank ECNs took 11.4% of global 
turnover. Inter-dealer trades, most of which are executed electronically, constituted 18.5% of turnover. The "dark 
pools" tend to be in the non-electronic execution segments of the market, with customer-direct transactions be- 
tween banks and end-user clients taking 24.4% of global turnover, and voice-broker trades accounting for 15.9% of 
global turnover (BIS, 2010, p. 16). 

There is also evidence of "home bias" in FX transactions. FX trades initiated in the US and Canada tend 
to take place during North American trading hours, those initiated in Japan and Australia tend to occur during 
Asian trading hours, and so on (D'Souza, 2008, Table 2). Estimation of impulse-response functions for the way 
exchange rates respond to local order flow data indicate that "dealers operating both at the same time and in 
the same geographic region as fundamentally driven customers have a natural informational advantage" (D'Souza, 
2008, pp. 23-24). In a similar vein, Covrig and Melvin (2002) found that Tokyo-based traders had an informational 
advantage over foreign-based traders regarding the course of the Japanese yen exchange rates. 

In previous papers, Kozyakin et al. (2010) and Cross et al. (2012), the implications of this fragmentation 
of information on FX markets were analysed by a combinatorial analysis of the different possible sequences of 
arbitrage operations. The "home bias" asymmetry in information is represented by having FX trader-arbitrageurs 
initially know only the exchange rates involving their own domestic currency. The question is then whether 
triangular arbitrage transactions are profitable because the cross-exchange rates are misaligned in relation to the 
exchange rates for the currency pairs. The arbitrage sequences that will be pursued would depend on which trader- 
arbitrageur first discovers the mis-alignment of the cross-exchange rates. The arbitrage operations are triangular, 
so, for example, if the US dollar trader were to discover that the euro-sterling rate was out of line with exchange 
rates for the dollar-euro and dollar-sterling currency pairs, an arbitrage profit could be made by selling dollars for 
sterling and using the sterling to buy euros. 

In the 3-currency case arbitrage operations are reasonable straightforward. The order in which trader-arbi- 
trageurs discover information about cross-exchange rate discrepancies makes a difference to the arbitrage sequences 
that will be pursued, and to the resulting new no-arbitrage ensemble of exchange rates that will emerge, but one 
arbitrage transaction suffices to eliminate the arbitrage opportunity. The 4-currency case is significantly more 
complicated, there being A P% = 24 possible triangular arbitrage operations. In Kozyakin et al. (2010) and Cross 
et al. (2012) it was shown that in this 4-currency world arbitrage sequences tend to be periodic in nature or 
display exponential behaviour, showing no tendency to approach a no-arbitrage ensemble of exchange rates in 
which there are no profitable opportunities. The intuitive explanation for this periodicity is that the exploitation 
of one arbitrage opportunity has "ripple effects," disturbing the ensemble of exchange rates and so creating further 
active arbitrage opportunities. 

The task of the present paper is to extend the analysis of FX arbitrage operations to a world of (^-currencies, 
where d is the number of currencies and trade-arbitrageurs involved. One key finding is that in the (d > 5)-currency 
case the surprising result is that the arbitrage sequences may follow a double exponential process as well as the 
periodicity or exponential behaviour observed in the 4-currency case. There is thus an "inheritance of instability" 
as we move to higher-order currency worlds. 

The structure of the paper is as follows. In Section 2 a mathematical formulation of the problem is given. In 
Section 3 it is shown that the arbitrage dynamics, initially specified in terms of some nonlinear operations, can be 
described by asynchronous matrix products. For the general case of an arbitrary number of trader-arbitrageurs the 
properties of the resulting matrices are studied in Section 4, The exposition of Sections 2 and 3 follows the work 
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of Kozyakin et al. (2010). In Section 5 basic results from Kozyakin et al. (2010) and Cross et al. (2012) about the 
arbitrage dynamics for an FX currency market with 4 traders are recalled. In particular, these results demonstrate 
that the exchange rates in a foreign exchange currency market with 4 traders, under appropriate choice of the 
arbitrage sequences, display periodic or exponential behaviour. In Section 6 we construct an example showing that 
in the case of 5 or more arbitrage traders the situation may be, in a sense, even worse — in this case the exchange 
rates may change not only periodically or grow exponentially, as in the case of 4 traders, but they may grow in 
accordance with a double exponential law. In Appendix A and Appendix B an explicit form for all the matrices 
involved in the description of the arbitrage dynamics for the cases d = 4 and d — 5, respectively, is presented. 



2. Statement of the Problem 

Consider a foreign exchange (FX) currency market involving d currencies involving the exchange rates of 
the i-th currency to the j-th currency, j ^ i. Naturally, the value of rji, the exchange rate of the j-th currency to 
the i-th currency, is reciprocal to the value of : 

1 

r ij 

The exchange rates vary with time depending on the state of the market and on the relations between the 
exchange rates. For instance, if at some moment the trader of the currency i realises that exchange rate of the z-th 
currency to the j-th currency, via the intermediate currency fc, may bring a profit, i.e., 

rik ■ r kj > r^, ij^j, i,j, 
then he can establish a new exchange rate for the currency 

So, in what follows, we will suppose that the exchange rates in our arbitrage system are updated in accordance 
with the following law: 

^ij.ncw max • Tfcj ■, ^jj ; old} ) 

1 

^ij,new 

and simultaneously the exchange rates for only one pair of currencies may be updated. 
By introducing the auxiliary quantities 

<Hj = logr-y, Vi ^ j, 

it is possible to pass from the "multiplicative" statement of the problem about arbitrage dynamics given above to 
the "additive" statement of the problem, which in this case will look as follows. Given a skew-symmetric d x d 
matrix A = (ay), for a triplet of pairwise distinct indices (i,j,k), i ^ j, k ^ i,j, the elements a.y and are 
updated in accordance with the following law: 

flij.ncw = maxjoifc + a>kj, aij j0 id} , (1) 

^'■i,ncw ^ij,new (^) 

The triplet of indices to = k), i ^ j, k ^ will be called the arbitrage rule. 



3. Linear Reformulation 

Since by ( 1 ) and ( 2 ) the element ct^ of the matrix A is updated "simultaneously" with the updating of the 
symmetric element aji, then it is reasonable to speak about updating of the pair (ay, a>ji)- Then, according to (1) 
and (2), the pair of elements (a^, aji) can be updated by one of two following scenarios. Either at first the element 
aij is updated in accordance with ( 1 ) , and then the element aji is updated in accordance with ( 2 ) , or at first the 
element aji is updated by the formula 

Oji,ncw = Uiax.{ajk + Ofci, dji. old} ; 
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which, in view of the skew- symmetry of the matrix A, leads to the formula 



a 



jz,ncw 



maxf-ay - a ik , - a ljMd } = - wm{a ik + a kj , Oij tM } , 



and then the element a^ is updated as follows 

Oij,ncw — %'i,riewi 

which amounts to the final expression for ciy. ncw : 

a>ij,new = min{aifc + ajy, ay,oid} ■ (3) 

Formulae (1) and (3) show that in the process of updating of the pair (fly, fflj-i) the element ay, independently 
of the scenario of updating, either is not changed (and then the matrix A is not changed, too) or this element may 
change its value as follows: 

^z?,new &ik ~\~ &kj : (4) 

whereas the symmetric element aji changes its value in accordance with (2) which is equivalent to 

Oji,ncw = a jk + a ki- (5) 

So, the following lemma is proved. 

Lemma 1. If the arbitrage law (1) and (2) holds then any pair of elements (aij,aji) either is not changed or is 
changed by the linear law (4), (5). 

Clearly, the converse assertion is also true. 

Lemma 2. If a pair of elements (a%j,aji) were to be updated by the law (4), (5), i.e., 

ya^j^CLjij I y (&zj',ncw j ^ji,ncw) 

then one of two scenarios of updating of the pair of elements {aij, aji) can be selected ( at first updating aij by formula 
(1) and then adjusting aji by the formula (2); or at first updating aji by the formula (1) and then adjusting a™ by 
the formula (2) ), under which the same pair of new elements (ay !now , aji jllew ) will be obtained. Namely, 

• if aij^ ncw > a^ after updating the pair {aij,aji) then, to interpret such an updating by the arbitrage law (1) 
and (2), one should first apply formula (1) and then formula (2); 

• i/oij.now < aij after updating the pair (aij, aji) then, to interpret such an updating by the arbitrage law (1) 
and (2), one should first apply formula (2) and then formula (1); 

• and last, ifaij^ ncw — Oy then the order of updating the elements of the pair (aij, aji), as well as whether such 
an updating took place at all, is inessential. 

The previous lemmata indicate that under the max-statement of the problem, as well as under its linear 
reformulation, it suffices to investigate only the dynamics of the pairs (aij,aji) of mutually symmetric elements 
since such an analysis allows us to reproduce also the scenarios of intra-pair updating of the elements a^ under the 
max-statement of the problem. 



4. Asynchronous Matrix Products 

Define in the space of all skew-symmetric d x d matrices A — (ay) a basis {e^}, 1 < i < j < d, by enumerating 
in some order the upper off-diagonal elements of such matrices. For example, let us define the element e^ , 
1 < i < j < d, of this basis as the skew-symmetric matrix whose element a^ is equal to 1 and the others vanish. 
Then each matrix in this basis can be represented as a column- vector: 

x = {<zi2, ai3, . . . , ai d , a.23, a,24, ■ ■ ■ a 2 d, ■ ■ ■ , a d -i,d} T € M d(d_1)/2 , (6) 
and computation of the new vector x ncw by the old one a; id will be defined by the expression 

*^ncw -^ujXqI^. (7) 
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Here the matrix B u is determined by the arbitrage rule u = k) applied to this step of updating and has the 
form typical in the theory of asynchronous systems, see, e.g., Asarin et al. (1992): all of its rows except one coincide 
with the rows of the identity matrix while one row contains exactly two non-zero elements specified by formula (4). 

If we define on the space of vectors (6) the usual Euclidean inner-product (•,•), and then use the matrix in 
(7) corresponding to the arbitrage rule w = (i,j,k), we will get the following representation 1 : 

{x + (-eij - e ki + e kj ,x)eij for k <i < j, 
x + (-e^ + e ik + e k j,x)eij for i <k < j, (8) 
x + (-eij + ejfc - ejk,x)eij for i < j < k. 

Since all the vectors (6) are column- vectors, then the relations (8) can also be rewritten in the following matrix 
form: 

f 1 - e y ( e fj + e L ~ e lj ) for k<i< j, 
B U =<I- e l3 (efj - ej k - e^) for i < k < j, (9) 
[ 1 - e y ( e I 3 ~ e Ik + e jk) for i<3 <k, 

where the upper index T denotes transposition of a vector or a matrix. 
Let us make some remarks resulting from formulae (8) and (9). 

Remark 1. Each matrix B u is a matrix with integer entries. This means that, when considering convergence 
issues for the matrix products involving such matrices, the convergence, provided that it takes place, will be always 
achieved by a finite number of steps. □ 

Remark 2. Each matrix B u is a projector. □ 

Remark 3. Since the procedure of seguential updating of vectors x according to (7) linearly depends on the initial 
vectors, then to analyse their convergence it suffices to consider only the initial vectors taking the values of the 
standard basis vectors {e, 3 -}. Then the whole problem becomes an integer-valued problem, i.e, a combinatorial 
one. □ 

Remark 4. The procedure of constructing the matrices B u is very similar to that of constructing the so-called 
"mixtures" of matrices under investigation in asynchronous systems Asarin et al. (1992). The difference is that in 
Asarin et al. (1992) each coordinate of the state vector x new in (7) is updated by a single rule, whereas in our case 
each coordinate of the state vector x new in (7) may be updated in accordance with several rules — see Section 5 
below. This is explained by the fact that the set of all the arbitrage rules •) having the same first two indices 
in general, contains more than one element. 

Investigation of convergence of the procedures (7) has very much in common with investigation of the joint/gen- 
eralized spectral radius of the family of all the matrices B^, see, e.g., Jungers (2009) and bibliography therein. 

So, a great variety of results from the theory of asynchronous systems and from the theory of the joint/generalized 
spectral radius might be helpful here. □ 

Clearly the arbitrage procedure (7) may have limiting fixed states (with respect to all the matrices B^) if and 
only if the eigenspaces of the matrices B u , corresponding to the eigenvalue 1, have nontrivial intersection F. To 
find this subspace F it is necessary to solve the following system of linear equations: 

x = B^x, V uj = (i,j,k) : 1 < i < j < d, k ^ 

Lemma 3. The subspace F of common fixed points of all the matrices {B^} consists of the column-vectors x = 
(xij) T satisfying 

Xij — —X\i + xij, 2 < i < j < d, 

where x\2, X\$, . . . , x\d are free variables. 

Proof. By induction with respect to d > 3. □ 

1 Recall that always i < j and k 
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Remark 5. By Lemma 3, one can assert that a skew- symmetric matrix A is unchangeable by application of any 
arbitrage rule if and only if it has the following form: 



( a 12 ai 3 ai4 

dl3 - <2i2 ClU - 012 

ou - a>i3 




V 



0>l,d-l did \ 

ai,d-i ^ a i2 aid — a i2 

a-\,d-i — «i3 o-id — ai3 

a l,d-l — a 14 Old — a 14 







aid ^ ai,d-i 




where the subdiagonal elements are defined by skew-symmetry and so are ignored. 
Lemma 4. The matrix 



p = E e ^ e ij+ E e a( e Tj- e u) 



□ 



(10) 



2<j<d 2<i<j<d 

is a projector on the subspace F. The vectors 

fi = Pen = ^ e J l ~ X! e ^-' * = 2, 3, . . . , d, 

l<j<i i<j<d 

where the second sum is assumed to vanish when i = d, form a basis ofV. 
Proof. Follows from Lemma 3 and Remark 5. 



□ 



Let us define in the space ]R d ( d a new system of coordinates by passing from the vectors x — (xij) T to the 

T 



vectors y = (yij) in accordance with the rule 



for 1 = i < j < d; 



Xij + Xn — Xij for 2 < i < j < d. 

Then the inverse change of variables is defined as 

_ \Vij for 1 = i < j < d; 

\ Uij ~ Vn + Vij for 2 < i < j < d. 



(11) 



(12) 



Remark, that in terms of the variables y, the subspace F of common fixed points of the matrices {B u } can be 
characterised as follows: 

F = {y = : Vij = for 2 < i < j < d} . 

Relations (11) and (12) can be rewritten in matrix form 

y = Q^x, x = Qy, 

where 



Q = I+ e ^( e y-e^), Q 1 =I- e n( e Ij ~ e u)- 

2<i<j<d 2<i<j<d 

Since the identity matrix / allows the representation 



(13) 



l<i<j<d 



T _ 
e ij e ij ~ 



E 

2<j<d 



T 



E ei i e S' 

2<i<j<d 



then by (10) we have 



q = e f v ; /; + p - 

2<i<j<d 



Therefore 

/ 4 = Pen 



)e u , i = 2,3, ...,d, 



and then 

e li = Q~ 1 f i , i = 2,3,...,d. 
Because each of the vectors i = 2,3, . . . , d, is a fixed point of each of the matrices B ull i.e. 

ft = B^fi, i = 2,3,...,d, Vw, 
then each of the matrices Q~~ 1 B U1 Q takes the block-triangular form: 



D u := Q- l BM = y Q £ J , (14) 
where 

en = -D^eii, i = 2, 3, . . . , d, V w, 

To describe the structure of the matrices D u in more detail, make use of the representations (9) and (13), and 
of the fact that e^e mn = 1 if and only if i = m and j = n, whereas e^e m „ = in other cases. 

Case 1 = k < i < j. By (9), in this case B u = I — eij(efj + - e^), from which 

D u = Q- X BM 

= (j — £mn( e ln ~ e lm)) (j ~ e ij( e Ij + e li — e lj") J (j + e mn(eJ n — e lm) 

2<m<n<d 2<m<n<d 

v ^ ~ e mn( e ln ~ e lm)) ( e ij( e ij + e li ~~ e lj)) + e mn( e l„ ^ e l ra ) 

2<m<r, 

/ — eij(efj + — ey) — ey(e^ + — e^) e mn (ej n — ef m ) 

2<m<n<d 

= I — e.ij(efj + e xi — ey) — e^e^- — e 14 ) = I — e.^e\ 



2<rn<n<d 2<m<n<d 



I e ij( e ij + e li e lj)(^^ e r?j.ri( e l n e l m , 

2<m<n<d 



Case 1 < k < i < j . By (9), in this case £? w = J — eij(efj + — ejL), from which 



A, = Q-^Q 

^ — ^2 e mn{eJ n — SlnSj (l — e ij( e fj + e ki ~ e kj)j (l + e mn( e ln — e l 

2<m<n<d 2<m<n<d 



2<m<n<d 2<m<n<rf 

J — eij(efj + — ejQf/ + ^ e m „(ef n — ef m ) 



2<m<n<(i 

J — eij(efj + e£i — ejQ — eij(ejj — e^) — e^e^ — e^.) + e^e^- — ej k ) 

= I — eij(e[j + eJi — <-/ , i 
Case 1 — i < k < j . By (9), in this case B^ = I — ey(ey — ef fc — ejL), from which 

2<m<n<d 2<m<n<d 

= I — (j — X e m „(ef„ — ef TO )J fe^^ey — ej k — e^j)) (j + ^ e mn (ef„ — ef TO ) 

2<m<n<d 2<m<n<d 

= I-(l- e mn(ei n - &f m )J (e y (e y - - e^) - ei^e^ - ef fe )J 

i 

= /+(/- X e mn (ef„ - eL) eye| = / + [ e xj + ^ e jn - ^ e mj ) e 



2<rri<n<d 



2<m<n<d j<n 2<m<j 



Case 1 < i < k < j. By (9), in this case B u = I — e i: ,(e^ — ef k — ejy), from which 



I e m „(e^ n e im)) [J e ij( e Ij e Ik e fcj))(^+ e„ m (ef n e 

2<m<n<d 2<m<n<d 

7 — ^7 — e m „(ef„ — ef m )^ feij(ey — ef k — ejQ^ + ^ e„ m (ef„ — e 

2<m<n<d 

7 — eijle^ — e ik — e fej ) ( 7 + e mn (e ln — e lm ) 



T • 
lm. 



2<m<n<d 

7 — eij{efj — ej k — ejy) — ey(e^ — e^) + ey(e lfc — e^) + e i:) (e^ — e^ fe ) 



' e ij( e Ij <'!>, e kj) — • 



Case 1 — i < j < k. By (9), in this case S w = 7 — eij(efj — ej k + from which 
7> w = Q-'B^Q 

)(j-e lj (e[ j -e[ k + ej k yj(l+ ^ e ron (ef n - ef m ) 

2<m<n<d 



-('- 


^ ] e mn( e l n e 




2<m<n<d 


= 7- 


7 — ^ ] e mn( e ln 




2<m<n<(i 


= 7- 


7 — ^ ^ e mn( e ln 




2<m<n<rf 







In c lm 



In c lm 



))( e y( e ij - e ffe + e Jfe)) ( J + e ™«( e i«- e i™) 

2 < m < n < d 

)) ( e ij( e ij - e ife + e Jfc) + eijOifc - e ij) 



2<m<n<d j<n 2<m<j 



T 

jk- 



Case 1 < i < j < k. By (9), in this case 7?^ = 7 — ey(e^ — ef k + ef k ), from which 



7 5^ e mn( e l„ e lr?j.)) (7 e ij( e ij e ik + e jfc)) M + e mn (e ln S.\ m , 

2<m<n<d 2<m<n<rf 

7 — ^7 — ^ e mn(ein — e i m )j ( e ij( e fj ~ e Ik + e Jfc)) M + ^ e m „(ef„ — ef m/ 

2<rn<n<d 2<m<n<rf 

7 — e ij( e fj — e fk + e Jfc) (7 + e mn(ef„ — ef m ) 



2<m<n<(i 

= 7 — eij(efj — ef k + ej k ) — ey(ey — e^) + ejj(e lfe — e^) — e^e^. — e^) 



7 e i j(ef : j e[ k + ej k ) — /«', 



It is convenient to unite the obtained relations into a single formula: 



7> 



'7- 


T 










for 


1 = k < i < j; 


7 - 


e »j ( e ij 


+ el 


" ^) 






for 


1 < k < i < j; 


/ + 

< 


( e u + 




6_yn — 


E2<m<j e mj^ 




for 


1 = i < k < j; 


7- 


e ij( e fj 


T 

~ e ik 


- efc) 


= 5 W 




for 


1 < i < k < j; 


7 - 


( e u + 






E2<m<j e "»J / 




for 


1 = i < j < k; 


>7- 


e ij( e fj 


e ifc 








for 


1 < i < j < k. 



(15) 



Now, to compute the matrix F^, it is necessary to retain in the right-hand part of (15) only the summands of 
the form e mn ej ? in which m — 1 and n,p,q > 1. Similarly, to compute the matrix G u it is necessary to retain in 
the right-hand part of (15) only the summands of the form e mn el in which m, n,p,q > 1. In doing so one should 



8 



keep in mind that I = Ei< m<n <d e ™« e mn- Tms results in 



and 



for 1 < k < i < j, 

for 1 = i < k < j, 

for 1 < i < k < j, 

-eijej k for 1 = i < j < k, 

for 1 < i < j < k 



j — p-p T . 

Ej<n 



E 



2<m<j j ^kj 



1 -^3( e l 3 ~ e ik~ e k 3 ) 

U-ey(e£-e£. + eT fc ) 



for 


1 = k < i < j, 


for 


1 < k < i < j, 


e kj for 


1 = i < k < j, 


for 


1 < i < k < j, 


e Jk for 


1 = i < j < k, 


for 


1 < i < j < k. 



(16) 



Naturally, in (16), in contrast to (15), the identity matrix I acts in the subspace spanned over the vectors e r , 
2 < to < n < d, and therefore it is representable as I = J2 



2<m<n<d ^rnn^rnn- 



Remark 6. In formulae (15) and (16) the sums corresponding to an empty set of indices are assumed to vanish. 



For instance, ^ 



j<n e jn 



if 3 



E 



2<m<j c ™? 



if 3 = 2. 



□ 



Remark 7. The idea of highlighting the matrices G u is to lower the dimension of the studied matrices since, 
according to (14), the behaviour of the products of matrices G u essentially determines the properties of the products 
of matrices D u and . □ 



5. Case d — 4 

The detailed analysis of the FX currency market with 4 trader-arbitragers, i.e., when d = 4, was fulfilled in 
Kozyakin et al. (2010) and Cross et al. (2012). Recall some related results. In this case the matrix A is of the form 



A = 



V 





-an 
-0.13 



a.12 


-au 



a 13 014 ^ 

0<23 Q.24 

a 34 

-a 34 ) 



In terms of the upper off-diagonal elements, the arbitrage rules can be written as follows: 



' (123) 


^12, new 


0,13 


- 023, 


other elements unchanged, 


(124) 


^12, new 


O14 


— A24, 


other elements unchanged, 


(132) 


^13, new 


a\2 


+ 023, 


other elements unchanged, 


(134) 


^13, new 


O14 


— «34, 


other elements unchanged, 


(142) 


a 14,ncw 


ai2 


+ a 2 4, 


other elements unchanged, 


(143) 


a 14,ncw 


O13 


+ a 34 , 


other elements unchanged, 


(231) 


Q23,ncw = 


-a\2 


+ ai 3 , 


other elements unchanged, 


(234) 


G23,ncw = 


0-2A 


— a 34, 


other elements unchanged, 


(241) 


a 24,ncw = 


-ai2 


+ ai4, 


other elements unchanged, 


(243) 


024,ncw = 


a-23 


+ 034, 


other elements unchanged, 


(341) 


^34. new 


-Ol3 


+ ai4, 


other elements unchanged, 


. (342) 


^34, new 


-023 


+ a 2 A, 


other elements unchanged. 



The subspace of fixed points F consists of all the vectors x = (xij) T , coordinates of which satisfy 



£13 -2^3, aria = xu - £24, £13 
^12 + 2:24, ari4 = X13 + X34, x 23 

-X12 + X14, X 2 4 — X 2 3 + x 34 , x 34 



2:12 + £23, X13 

-X 12 +X 13 , X 2 3 

-X13 + XU, X 34 
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X14, - S34, 
2*24 - 2:34, 
-X23 + X 2 4- 








This system has redundant equations. In particular, all the equations containing the minus sign are redundant. 
By removing them we obtain: 

^13=^12+^23, Zl4 = X 12 + X 24 , ^14 = ^13+^34, X24 =X 2 3+X 34 . 

Here one of the equations determining x 44 is also redundant. By removing it, we obtain the final set of independent 
equations defining the subspace F: 

£13 = X12 + £23, xu = %13 + £34, 2:24 = £23 + 2:34 • 
By choosing X12, £13 and x 44 as independent variables, the remaining variables can be defined by the equalities 

2-23 = -Xl2 + 2^13, 2-24 = -X 12 + X U , X 34 = ~X 13 + X 14 . 

The form of the matrices B w , and G w for this case is given in Appendix A 

As was noted in Kozyakin et al. (2010) and Cross et al. (2012), in the case d = 4 the arbitrage process may be 
non-convergent for some periodical sequence of the arbitrage rules. To support this claim let us consider the set 
of all vectors each of which is mapped by some of the matrices to zero. Direct calculations show that, up to 
non-negative factors, this set consists of the following 12 vectors: 

si = I 1 I , s 2 = I I , s 3 = I 1 I , s 4 = I I , s 5 = I 1 I , s 6 

Sr = -Si, S 8 — -82, S 9 = -S 3 , S W = -S 4 , Sn = —S 5 , Sl 2 = —s 6 , 

Moreover, as is easily verified, the set 

§ = C0{±Si,±S 2 ,±S3,±S4,±S5,±S6} 

is mapped by each of the matrices into itself and is a body, i.e., contains interior points. Hence, S is a unit ball 
of some norm || • ||* for which 

||GJ* < 1, Vw. (17) 

So, the set of matrices {G^} is neutrally stable and, it seems, there is the possibility that all products of 
matrices G u are convergent. In Sections 5.1 and 5.2 it is shown that this is not the case. 

Remark 8. By (14) the following representation is valid 

B " = Q ( cl ) 

From here, and from the estimate (17), it then follows that the norm of any products of the matrices B u may 
grow no faster than cn with some constant c, where n is the number of factors B u in the product. Example 5.2 
demonstrates that such a rate of growth of products of the matrices B^, in the case d — 4, is really achievable. □ 

5.1. Example 

Consider the following sequence of the arbitrage rules 



(3,4,1), 


w(4) 


= (1,4,2), 


(2,3,1), 


w(8) 


= (3,4,2), 


(3,4,2), 


w(12) 


= (1,4,3), 


(1,2,4), 


w(16) 


= (1,4,3), 



and extend it by periodicity. Set also Xq — s%, and build recurrently the sequence 
x n = G u) ( n )X n ~i, n= 1,2, ...,16. 
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Then 



Xn 


= .Si 








X 1 


= G,, 


1 \ Xn = 


(t/i jo 1 ! = 


.So 


Xo 


= G,, 




( T ! 1 OQ'tiKl 

v - ' y\_Ao ) 1 




X'l 


= G, 




C T ( 'i A 1 \ — 


, 






4) x 3 — 


£?(142)3?3 = 


eg, 


f r 
A 5 




5) J -4 — 


Lt (134) j '4 — 


ft 4, 


k 6 


— G w 


6)^5 = 


^(243)^5 = 


«1, 


aj 7 


= G w 


7)^6 = 


G(231) a; 6 = 


-S5, 


a; 8 


— G w 


8)^7 = 


0(342)^7 = 


«3, 


CCg 


= G w 


9)^8 = 


G(241)^8 = 


S6, 


scio 


= G w 


10) #9 = 


= G(i34):E9 = 


= S4, 




= G w 


11)^10 


= G( 3 42)a;io 


= s 2 , 


aJi2 


= G u 


12)^11 


= G(U3) X 11 


= Si, 


^13 


= G w 


13)^12 


= G( 23 4)a;i2 


= si, 




= G w 


14)^13 


= G(i32)a;i3 


= «3, 




= G w 


15)^14 


= G( 124 )Xi4 


= -s 2 , 


ai6 


= G w 


16) #15 


= G(i4 3 )a;i5 


= -si. 













Continuing building, we will obtain a 32-periodic sequence x„. It is not difficult to see that the sequence of 
products of the corresponding matrices B^: 

H n = £>w(ri)-^ri-l, H = I, n = 1, 2, . . . , 

will be also 32-periodic and non-convergent to any limit. 

5.2. Example 

Consider the following sequence of the arbitrage rules 

w(l) = (1,4,3), w(2) = (3,4,1), w(3) = (3,4,2), w(4) = (1,4,2), 
w(5) = (1,2,4), «(6) = (2,3,1), w(7) = (1,3,2), w(8) = (2,4,3), 
w(9) = (1,3,4), w(10) = (2,4,1), w(ll) = (1,2,3), w(12) = (2,3,4), 

and extend it by periodicity. Set also Xq — si, and build recurrently the sequence 





— G^^^a;,!—! 


n=l,2, 


..,12. 










Xq 


= si, 






X\ 


= G w (i)cc = 


G(i43)a; = 


Sl, 


X 2 


= G u (2)a;i = 


0(341)03! = 


Sl, 


x 3 


= G w ( 3 )aj 2 = 


G(342)^2 = 


Sl, 


CC4 


= G w ( 4 )a; 3 = 


G(i42)3;3 = 


s 2 . 


x-, 


= G w ( 5 )a;4 = 


G(124)^4 = 


s 2 , 


cc 6 


= G^( 6 )a;5 = 


G(231)^5 = 


-s 6 , 


x 7 


= G w ( 7 )a; 6 = 


G(i32)a; 6 = 


-Se, 




= G u ^x 7 = 


0(243)^7 = 


-S3, 


X<) 


= G w ( 9 )a; 8 = 


G(i34):E8 = 


-si, 


3510 


— G w (io)^9 = 


= 0(241)^9 = 


= -S4, 


Xu 


= G w (n)a; 1 o 


= G(i23)a;io 


= -s 5 


3J12 


= G w ( 12 )a;ix 


= G(234)a;ii 


= Si. 











Continuing building, we will obtain a 12-periodic sequence x n . The sequence of the corresponding products of 
the matrices G^On) is bounded by (17) whereas the sequence of products of the matrices B u : 

H n = £ w ( n )if n _i, H = I, n = 1, 2, . . . , 

is divergent, and the norms of matrices H n have the linear rate of growth, see Remark 8. 

6. Case d — 5 

The full set of the matrices B u , D u and G w for this case is presented in Appendix B . Direct calculations show 
that the matrix 

H = £( 14 3)B(23i)B(245)-B(342)-S(451)-E?(124)-E?(453) 

is as follows 

/ 1 
1 
1 
000 

TT _ 1-1 

11 1 







\ 0-1 









-1 











\ 


























-1 


1 














1 


























1 














1 








1 




















1 














-1 


1 





























1 





1 

















/ 
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and its spectral radius p(H) is strictly greater than 1: 
P(H) = * 1.618. 

To ascertain this it suffices to note that the matrix H has the following eigenvalues (ordered by decreasing absolute 
values): 

1 + \/5 \/5 — 1 

i, i, i, i, i, , o, o, o. 

Consider the following sequence of the arbitrage rules 

u,(l) = (4,5,3), u,(2) = (1,2,4), 0,(3) = (4, 5, 1), w(4) = (3, 4, 2), 
o,(5) = (2,4,5), 0,(6) = (2, 3,1), 0,(7) = (1,4,3), 

and extend it by periodicity. Then, as is known from the linear algebra, for almost all initial vectors 2 xq, the 
sequence of vectors 



B u ( n )X n -i, n = l,2, 



obtained by sequential application of the arbitrage rules o,(n) will diverge with the rate I 1+ 2 j , i.e 



\ n/7 / \ n/7 

^ ||xo||<||*J|<C(i±^] ll^ol 



where c, C are some positive constants. Taking into account that 1.071 < ( ) < 1.072, the inequalities 



1/7 

o i 

2 

obtained may be represented in a more descriptive form: 

c 1.071" || x 1| < ||x n || < C1.072"||iCo||- (18) 

Remark 9. Relations (18) imply that for d = 5, m contrast to the case d = 4, ifte norms of the vectors {x n } as 
well as the norms of the matrix products 

H n — Bw(n)Hn-ii H — I, n = 1, 2, . . . , 

may grow with the exponential rate. □ 

Remark 10. As was shown in Remark 9, in the case d = 5, the norms of the matrix products may grow with the 
exponential rate whereas in the case d = 4, according to Section 5, the norms of the matrix products may have at 
most a linear rate of growth. Any seguence of arbitrage rules involving 4 traders may be treated as a particular 
case of a sequence of arbitrage rules for 5 traders. Therefore in the case d = 5 one may observe also the sequences 
of matrices whose norms of the products have a linear rate of growth, as well as sequences of matrices (vectors) 
varying periodically, see Section 5. 

When d > 5 we face an effect of "inheritance of instability" taking place in the cases d = 4, 5. In this case one 
may observe all the types of behaviour of the matrix products described above as well as of the corresponding state 
vectors x n : periodicity, growth at a linear rate and growth at an exponential rate. □ 

Remark 11. If we return back from the "artificial" quantities aij to the true exchange rates 

n 3 =e a *\ (19) 

we arrive to a rather unexpected conclusion. 

In the case d = 4 the exchange rates may vary periodically or grow exponentially which follows from (19) and 
the results of Section 5. 

2 To be more specific, for all initial vectors not belonging to the linear invariant subspace of the matrix H corresponding to the set 
of eigenvalues distinct from — ■ 
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In the case d > 5 the exchange rates also may vary periodically or grow exponentially but also, what is surprising, 
they may grow in accordance with the double exponential law. For instance, in the case d = 5 formula 
(18) implies 

e £1 - 071 " < K|| <e Gw72 \ 
where r n is the exchange rate vector 

{ru, r-ia, • ■ • , r ld , r 23 ,r 24 , . • . r 2c j, . . . , r d _i jd } T 

at the moment n. □ 

Remark 12. As follows from the analysis undertaken, after introducing new variables bringing the matrices 
to block-triangular form, one obtains a set of variables, which is updated only by matrices G u , and the behaviour of 
which is thus determined only by its own history. These variables may be called, in a sense, "leading exchange 
rates indices. " For d = A, this set is constituted of the variables 

fflia - ai3+a23, au - au + a 2 4 , ai3-ai4 + a 34 , (20) 
which, in terms of exchange rates, corresponds to the set of variables 
ri 2 -r 23 J"i2 • r 24 n 3 ■ r 34 

= r 12 ■ r 23 • r 3 i, = r i2 • r 24 • r 4i , = r i3 • r 34 • r 44 . (21) 

fi 3 r u r 14 

-For d = 5, i/ie related set is constituted of the variables 

a 12~ a l 3 + a 2 3 , a 12~ a l 4 + a 247 a 12 — a 15 + a 25j a l 3 ~ a 14 + a 3 4j a l 3 ~ a 15 + a 35: a 14~ a 15 + a 45i (22) 

and to it, in terms of exchange rates, corresponds the set of variables 

ri 2 • r 23 • r 3 i, ri 2 -r 24 -r 4 i, r 12 ■ r 25 • rsi, n.3 • r 34 • r 4 i , ri 3 • r 35 • r 5 i, ri 4 -r 45 -r 5 i. (23) 

Fbr d = A, the variables (20) and (21) always vary boundedly over time, while, for d > 5, t/ie variables (22) 
and (23) may (/row unboundedly. At the same time, the original set of exchange rates {r^ } and the corresponding 
set of "artificial" log-variables {a^ } for d — A, as well as for d > 5, may grow unboundedly. 

The intuition underlying the "leading exchange rates indices" in (23) can be explained as follows. The term 
ri 2 r 23 r 3 i can be interpreted in the following way: the 1st trader is selling the currency 1 for the currency 2, 
then uses the currency 2 to buy currency 3 and then uses the currency 3 to buy his "home currency" 1. So, the 
term ri 2 r 23 r 31 is the "arbitrage profit factor" for the 1st trader when the trades involve the "most simple" cyclic 
triangular path: 

H2->3^1. 

If this factor r 12 r 23 r 31 is equals 1 then such "cyclic triangular trades" bring no profit to the 1st trader, while 
when ri 2 r 23 r 3 i =/= 1 then such "cyclic triangular trades" bring profit (or loss) to the 1st trader. Hence the condition 

ri 2 r 23 r 3 i = 1 

is the "no-arbitrage" condition for trading by the cyclic path 1 — ► 2 — > 3 — > 1, for the 1st trader. This "unit 
circle" no-arbitrage condition was pointed out in Cournot (1929, first published in 1838). 

Then ( 23 ) is simply the full set of factors whose simultaneous equality to 1 represents the no-arbitrage conditions 
(NACs) for the whole set of currencies in the d = 5 case: 

ri 2 r 23 r 34 = 1, r i2 r 24 r 4 i = 1, ri 2 r 25 r 5 i = 1, r i3 r 34 r 4 i = 1, ri 3 r 35 r 5 i = 1, ri 4 r 45 r 5 i = 1. (24) 

It appears that different currencies are represented in (24) non-symmetrically. But in fact there are different 
NACs which are equivalent to each other. 
We may introduce a measure, say, 

Q = |ri 2 r 23 r 3 i - 1| + |ri 2 r 24 r 4 i - 1| + |ri 2 r 25 r 5 i - 1| + |ri 3 r 34 r 4 i - 1| + |ri 3 r 35 r 5 i - 1| + |ri 4 r 45 r 5 i - 1|, 

which can be interpreted as an "accumulative" NAC. The meaning is that the equality Q = involves the absence 
of profitable arbitrage opportunities in the whole FX market system, while the quantity Q ^ characterises the 
extent of profitable arbitrage opportunities in the FX market system. 
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7. Concluding Remarks 

The EMH implies that arbitrage operations in markets such as those for FX should be self-extinguishing. This 
simplistic view of arbitrage processes may reflect a failure to consider the complex chains of arbitrage possible in 
financial markets, such as those for FX, that have segmented trading platforms and fragmented information. In 
Kozyakin et al. (2010) and Cross et al. (2012) it was shown, applying combinatorial analysis and asynchronous 
systems theory, that arbitrage sequences in a 4-currency world tend to display periodicity or exponential growth 
rather than being self-extinguishing. The present paper has used a log-linear reformulation of the analysis to 
show that in a (d > 5)-currency world, as well as inheriting the periodicity or exponential growth found in the 
4-currency case, arbitrage sequences can be augmented by a double-exponential process. Thus there appears to be 
a cumulative "inheritance of instability" . 

The results may help to explain why extensive arbitrage trading is endemic in FX markets. The fact that the 
number of currencies makes a difference is interesting: periodicity or exponential behaviour arise when you move 
from 3 to 4 currencies; double exponential behaviour as well, when you move beyond 4 currencies. The Financial 
Times provides daily quotes for 52 currencies against the US dollar, euro and pound sterling. But some currencies 
are more important than others. BIS data indicate that trades involving the US dollar accounted for 85/200 of all 
FX trades in 2010, followed by the euro (39/200), Japanese yen (19/200) and sterling (13/200), so trades involving 
these four key currencies accounted for 78% of global FX turnover (BIS, 2010, p. 12). In terms of currency pairs 
the leading trades are the dollar-euro, which accounted for 28% of global turnover, followed by the dollar-yen 
(14%), dollar-sterling (9%), euro-yen (3%), euro-sterling (3%) and sterling-yen (3%) — see BIS (2010, p. 15). In 
the terminology of Section 6 of this paper, these four currencies may be called "leading currencies." In a world 
containing just these four leading currencies, arbitrage operations can display periodicity, thus being bounded, or 
"mild" exponential growth. Once we move to account for triangular arbitrage outside these four leading currencies 
the arbitrage sequences may be unbounded, displaying "explosive" double exponential behaviour. In practice 
central banks often intervene to attempt to prevent exchange rates exceeding certain bounds. The history of 
financial crises, however, is replete with examples of central bank interventions in FX markets being unsuccessful, 
the central banks being eventually obliged to abandon their exchange rate pegs or target ranges. 

The "leading exchange rate indices," discussed at the end of Section 6 of this paper, indicate arbitrage profit 
opportunities expressed in an FX trader's "home currency" unit of account. The August 15 1971 abandonment 
of the convertibility of the US dollar into gold at a fixed price of 35$ per ounce, and the breakdown of the 
Bretton Woods system whereby IMF member currencies were pegged to the US dollar, left the international 
monetary system without a well-defined numeraire or unit of account. To fill this gap, currency-specific effective 
exchange rate (EER) indices are calculated by national central banks and finance ministries, and by international 
organisations such as the IMF and BIS. An EER is a geometrically weighted average of the bilateral exchange rates 
for a particular currency. So, for example, the BIS currently publishes broad EERs covering 61 currencies, and 
narrow EERs covering 27 currencies, using 2008-2010 international trade flow data for weights (Klau and Fung, 
2006; BIS, 2012). The various EERs use the "home currency" unit of account to provide exchange rate indices 
designed for a world in which there is no gold or single reserve currency anchor. 

Special Drawing Rights (SDRs) were created by the IMF in 1969 to supplement the gold and foreign exchange 
reserves available to member central banks. The SDR was initially defined as equivalent to 0.888671 grams of fine 
gold, which could be exchanged for 1 US dollar in 1969, but is now defined as a weighted value of the US dollar, 
euro, pound sterling and Japanese yen, with the weights being changed every five years, expressed in US dollars 
(IMF, 2012). A European Currency Unit (ECU), conceived by European Economic Community (EEC, later the 
EU) in 1979, served as an EU unit of account, being a weighted value of EU member currencies, until the euro was 
launched in 1999. 

The governor of the People's Bank of China recently advocated that the SDR be adopted as a global, tradeable 
reserve currency (Xiaochuan, 2009). The SDR is a dollar-euro-sterling-yen currency index, expressed in US 
dollars, so this is akin to the "leading exchange rate" index faced by a dollar FX trader-arbitrageur in the d = 4 
case considered earlier in the present paper. Accordingly, the periodicity and exponential growth features of 
arbitrage sequences involving these four currencies would apply. This suggests that central banks holding the SDR 
as a tradeable reserve currency would be faced with endemic, rather than self-extinguishing, opportunities for 
arbitrage profits on their FX reserve portfolios. The proposal that "the basket of currencies forming the basis for 
SDR valuation should be expanded to include currencies of all major economics" (Xiaochuan, 2009, p. 3) would 
exacerbate the arbitrage problem, introducing the "inheritance of instability" arising in the d > 5 case considered 
in the present paper. 
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Appendix 



Appendix A. Matrices B^jk),D{ijk) an d G(ijk) in the Case d — 4 

The matrices B^: 
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The matrices Q, Q 1 and D u — Q 1 B UJ Q take the form: 
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Appendix B. Matrices -B(ijfc), -D(ijfc) and G^k) in the case d — 5 
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